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Abstract 

We study string-gas cosmology in dilaton gravity inspired by the fact that it naturally arises 
in a string theory context. Our main interest is the thermodynamical treatment of the string- 
gas and the resulting implications for the cosmology. Within an adiabatic approximation, 
thermodynamical equilibrium and a small, toroidal universe as initial conditions, we numer- 
ically solve the corresponding equations of motions in two different regimes describing the 
string-gas thermodynamics: (i) the Hagedorn regime, with a single scale factor, and (ii) an 
almost-radiation dominated regime, which includes the leading corrections due to the lightest 
Kaluza Klein and winding modes, with two scale factors. The scale factor in the Hagedorn 
regime exhibits very slow time evolution with nearly constant energy and negligible pressure. 
By contrast, in case (ii) we find interesting cosmological solutions where the large dimensions 
continue to expand and the small ones are kept undetectably small. 



1 Introduction 



String cosmology is rapidly growing in importance for at least two reasons. First, unless there is 
a fortunate conspiracy of scales and one or more of the extra dimensions are ~ 0(Tev _1 ) PP, we 
cannot reasonably expect to detect any stringy effects in accelerators in the foreseable future. 
Hence string cosmology may provide the only means of testing string or M-theory concretely, 
through e.g. non-commutative effects 2 a or the AdS/CFT correspondence [Sj. 

On the other hand, traditional cosmology has been fettered by the chains of the singularity 
theorems of General Relativity and is therefore usually only seen as understanding the evolution 
of our universe over the past 15 billion years or so since the Big Bang. String theory offers the 
exciting possibility of a resolution to the big bang singularity thereby opening up a potentially 
infinite pre-history of the big-bang. Models constructed in this vein include the pre-big-bang 
[UE] (see also H3IZ!) and the recent ekpyrotic/cyclic models [HI El which has lead to new work 
on string propagation in orbifold backgrounds with curvature singularities |10| . 

Traditional cosmology, because of the limitations of General Relativity, suffers from another 
great lacking: it is unable to make any predictions about the number of dimensions we live 
in or about the spatial topology of the Universe. String and M-theory theory, in contrast, 
predict that we live in either 10 or 11 spacetime dimensions. Perhaps the greatest challenge for 
string-cosmology, after understanding the big-bang singularity, is to explain why and how three 
dimensions became observable and large while 6 or 7 are either small or unobservable for some 
other reason. Perhaps the only proposal in this direction so far is the so-called Brandenberger- 
Vafa (BV) scenario jllj where it is assumed that the universe is small and compact and that 
exactly three space dimensions become large because of the dynamics of winding modes, which 
play a particularly important role. Subsequently it has been pointed out in [I] that the low 
energy effective action of a dilaton-gravity system, naturally emerging in string theory, has a 
duality symmetry that is a manifestation of the string T-duality R — > a' /R symmetry, that 
plays a crucial role in the analysis of |llj . In this respect, the dilaton-gravity system is more 
suitable than standard General Relativity for the BV scenario Although some works have 
already shown that the BV scenario can be realized and extended in a variety of ways |12| I13| 
El El El El El El 1201 1^1 E2] (see also j2B] ) , a full string-theory analysis is very complicated 
and still not completed, even in the simplest toroidal compactification. 

In this paper, inspired by the BV scenario, we numerically solve the dilaton-gravity equations 
of motion ^2] with some types of "stringy matter" . Adiabatic evolution (which implies constant 
entropy), weak string coupling and thermodynamical equilibrium are always assumed in our 
analysis. For simplicity, we analyze Type IIA or IIB closed string theory on a T 9 torus, with no 
branes 1 . In particular, we consider the following two regimes: 

• (i) Hagedorn matter at high energy densities in a very small homogeneous and isotropic 



universe with a common compactification radius ~Va', and 



• (ii) an almost-radiation dominated regime with two independent scale factors, associated 
with the large and small dimensions. 

In the latter case, the lightest Kaluza Klein (KK) and winding mode contributions are also taken 
into account. 

1 In |17| it has been shown that fundamental strings, even in the presence of D-branes, are still the dominant 
degrees of freedom for the realization of the BV scenario. 
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It is important to stress here that in both of these regimes, our matter is manifestly T- 
duality invariant. This symmetry is broken in our set-up only by the (arbitrary) choice of initial 
values. The main relevant questions in the two cases are respectively: which is the evolution 
of the universe at early times for a free Hagedorn string gas in thermal equilibrium? Assuming 
large and small dimensions as initial conditions, how do they evolve? In particular, do the small 
dimensions remain small? 

String matter in the Hagedorn phase was already briefly discussed in ^2] where it was realized 
that, to a first approximation in which the energy is constant, it leads to a very slow evolution of 
the universe, as in eq. (|24|) . We complete the analysis of ^2] by relaxing this approximation and 
imposing the conservation laws for KK and winding modes. We find that in any practical sense 
there is no departure from the behaviour dictated by eq. (|24[) . and thus no relevant dynamics 
emerges in this set-up. It is important to recall how to interpret this result and its connection 
to the original BV scenario. The above result is obtained by assuming thermal equilibrium and 
a free, ideal, string gas, whereas the dynamics and interactions of winding modes at very early 
times are crucial in the BV scenario. The results we get in case (i) are therefore not closely 
connected with the original BV proposal 11] or, rather, they have to be taken the other way 
around. Namely, unless string interactions are taken into account and/or thermal equilibrium 
is relaxed, no interesting dynamics emerges. 

On the other hand, the evolution of the small and large dimensions in case (ii) is much 
more interesting. First of all, we will show that when there is only pure radiation, the small 
dimensions can be stabilized and kept small relative to the large dimensions. Essentially, it is 
only required that the initial expansion rate of the large dimensions is bigger or of the same 
order than that of the small ones 2 . This mainly comes from the fact that the pressure in the 
small dimensions vanishes in the case of pure radiation. 

When matter, in the form of KK and winding modes, is included, the choice of initial 
conditions becomes more relevant. The crucial point is played by winding modes that are able to 
distinguish large and small dimensions, leading respectively to a positive/negative contribution 
to the pressure along the large/small dimensions, as is clear from eqs. (|46|) and (|47|) . This turns 
out in an expansion of the large dimensions and at the same time the possibility of keeping 
almost constant the small ones. In fact we found that there exists a wide range of parameters 
for which the small dimensions actually remain small (see fig. 0} , while the large ones expand 
as required in the presence of radiation and string matter (see fig. [SJ). This is actually achieved 
for the natural initial condition where the small dimensions are close to the the self-dual radius 



V" , together with the condition of the expansion rate mentioned in the pure radiation case. 



These are the most important results of the paper. There exist a wide range of initial con- 
ditions for which the small dimensions are stabilized around the self-dual radius before entering 
into a purely radiation dominated phase, regime in which they are asymptotically stabilized 
anyway to a nearly constant value. 

There are clearly several open issues that we do not consider but that deserve further in- 
vestigation. First of all, thermal equilibrium, as well as isotropy and homogeneity along all 
directions, should be relaxed and string interactions need to be taken into account. Further, 
since toroidal compactifications are not realistic, always leading to non-chiral theories, study of 
more general topologies is required. See |191 12U| for recent progress in this direction. 

2 It is actually quite hard to imagine initial conditions where the expansion rate of the small dimensions are 
bigger than that of the large ones. 
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Another very important issue unexplored in this work is the transition from dilaton gravity 
to General Relativity. A common lore is to assume that supersymmetry breaking occurs at later 
times and stabilizes all moduli, including the dilaton, which then becomes massive 3 . It is not 
clear how the evolution of the large and small dimensions are altered by this transition to an 
effective description in terms of General Relativity. Notice, however, that if one considers in case 
(ii) pure radiation with no matter, then the late-time behaviour of the dilatonic gravity solutions 
lead to a standard Friedmann-Robertson- Walker (FRW) radiation-dominated cosmology for the 
large dimensions, whereas the small ones are constant. This behaviour nicely matches with the 
solution found by 20 in General Relativity in a purely KK extra-dimensional scenario. 

Finally, we assume an adiabatic approximation to simplify the analysis. A deeper justification 
of this approximation, or the dynamics which arises when it is dropped, is required. One can 
naively expect dropping this condition to yield a wider variety of dynamics, including more 
violent evolution. 

This paper is organized as follows. In section 2 some general aspects of dilaton gravity are 
briefly reviewed. In section 3 we analyse the dynamics of the system in the extreme Hagedorn 
regime of high energy densities with a single scale factor. In section 4 the dilaton-gravity 
equations are solved in an almost-radiation dominated regime with large and small scale factors, 
as well as in the presence of some massive stringy matter. We present some conclusions in the 
final section and present some useful formulae in the appendix. We set the string scale a' = 1 
in the following. 



2 Dilaton gravity 

We shall study the dilaton gravity equations of motion with a massless dilaton field corre- 
sponding to the low-energy effective action of string theory in D + 1 space-time dimensions, 
described by [103H2]: 

e - 2<s >{R + 4(V<i>) 2 }+£ M ], (1) 

where g is the determinant of the background metric gw, and Cm corresponds to the Lagrangian 
of some matter. The coupling of with gravity is the standard one arising in string theory. 
Hereafter we shall consider the case with D = 9. We are interested in the case where the whole 
universe is small and compact, corresponding to a flat T 9 -torus. In this case, if one considers 
field configurations that are spatially homogeneous, assuming only slow time-dependence and 
adopting an adiabatic approximation it has been shown in E] that the action Jy) exhibits a 
duality symmetry, a low energy manifestation of the string T-duality R — > 1/R symmetry. The 
ansatz for the metric and dilaton we use is 

9 

ds 2 = -dt 2 + R 2 i{t)dx 2 , Ri = e K{t) , $ = $(t) . (2) 

i=l 

Here the Ri denotes the i-th scale-factor of the torus. The equations of motion simplify if one 
introduces a shifted dilaton, ip, via 

9 

il> = 2$ - £ A, . (3) 

i=l 

3 An alternative would be the models in which the dilaton stays massless but has universal couplings. It is even 
possible to reconcile this with quintessence if the dilaton has a divergent runaway behaviour |24l I25| . 



S = / d D+1 x^g- 
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Given the metric © the equations of motions of the dilaton-gravity system are then |12j 



-EA? + V> 2 = e+E, (4) 
i=l 

\i-i>\i = \e^P h (5) 

^-EA t 2 = \e+E, (6) 
i=i / 

with E the total energy and Pi the total pressure along the i-th direction found by multiplying 
the total spatial volume of the space by the energy density and pressure appearing in Cm of (JJ). 
Here a dot denotes derivative with respect to cosmic time, t. These equations are manifestly 
invariant under the duality symmetry [3J El 

Aj — > — Aj, *-*-£Ai, (7) 

i 

under which -0 defined in eq. (j3J) is left invariant. It is typically assumed that the scale factors Ri 
are the same in all directions, i.e. Ri = R. In contrast, we also consider in section 4 a scenario 
where the background is homogeneous and isotropic in cf-spatial large dimensions and (9 — d)- 
spatial small dimensions. We denote the large and small dimensions with their corresponding 
scale factors, as 

R = , r = e u . (8) 
In this case eqs. Q-© take the form 

- dfi 2 - (9 - d)v 2 + ^ = e*E, (9) 
fl-iPfi = ±e*P d , (10) 

i>-fc = \^Ps-d, (ii) 
4> - dfi 2 - (9 - d)v 2 = -e*E, (12) 



where 



dF 

P d = - — , Vi = l,...,d, (13) 



dm 

P^d = ~, Vi = d+1,...,9, (14) 



in terms of the free energy F. 



3 Hagedorn regime 

According to the original BV proposal the very early universe was compact over all nine 
spatial dimensions with radii r ~ 1 in string units. In this section we study the dilaton- 
gravity equations with E and P the total energy and pressure of a free string gas in 
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thermal equilibrium in such a compact, small universe. For simplicity, we consider the string gas 
associated with type IIA/IIB string theory compactified on a square T 9 -torus, simple product of 
nine circles, with radii all equal to a common value r = e x . Although the strict thermodynamical 
limit V — > oo cannot be taken for this system, thermodynamics is still trustable as long as the 
system contains many degrees of freedom. In our case, this implies having an energy density 
p S> 1. Our first step is then to derive the equation of state of the string matter in this regime, 
or equivalently the energy and pressure entering in eqs. (HJ)-©. 



3.1 Microcanonical ensemble 



String thermodynamics in the high density phase has been a subject of intense study in the past 
(for a review see e.g. |27|). The most important point is the emergence of a critical temperature, 
the Hagedorn temperature Th |28| . where the partition function of a free string gas diverges. It 
was soon realized that in this regime the usual thermodynamical equivalence between the canon- 
ical and microcanonical ensembles can break down and the latter, more fundamental ensemble, 
must hence be used. Furthermore, due to the presence of winding modes, the thermodynamical 
properties of a string system at finite volume differs substantially from that of an infinitely 
extended system ^2 EH]- Following j22j, it is useful to derive the energy density of states Q(E) 
by analyzing the singularities of the one-loop string partition function in the complex /3-plane, 
where (3 = 1/T is the inverse temperature. This analysis has already been carried out in [22] 
but for completeness we report some useful details in the Appendix. 

For a totally compact space at high energy, the leading singularity (a simple pole) of the 
partition function at (3 = [3h is not sufficient to establish the thermodynamical properties of the 
system 129j . The first next-to-leading singularities are poles of order 18 with a dependence 
on the compactification radius. In this case a useful parametrization of the one- loop partition 
function Z is given in eq. (|57|). from which one computes the density of states Q(E) by means 
of eq. (|58|) and thus the associated entropy S = \ogQ(E) 4 . It reads 



S(E,r)~p H E + logll 



1 



r(is)^V 



{n K E) 17 rj$e-™ E + ( m E) x \fe-^ E 



(15) 



where /3jj = 2\/27r, and 




(16) 



2 - ^4 - 2r 2 



The energy as a function of r is given directly by (|15j) . since S = constant, by the assumption 
of adiabatic evolution. On the other hand, the temperature and pressure, defined as 



1 

T 



dS_ 

BE' 



OS 



9 <9(logr) 



(17) 



4 We have numerically checked that the term A(/3, R) in eq. 15711 is negligible and thus is not reported in the 
following. 
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yield 



1 

T 



1 



1 



where 



T H xET(18) 

Til 

9x\' 



1 



18 



y 17 (17 -y)e~y + 



is 



z 17 (17-z) 



z e 




(18) 
(19) 

(20) 



and y = rjxE, z = rjwE, w = t]kwE. As usual in the microcanonical ensemble, the temperature 
is a derived quantity (from S and E) and its explicit form is needed only to compute the pressure 
P. When the radius r is close to unity, T and P are approximately given by 



l.^ 7 L + Cl E 17 e 
T T H 



-ijE 



P 



C 2 E 17 e^ E 



(21) 
(22) 



where r) ~ rjx ~ r/w for r ~ 1, and C\ and C 2 are certain polynomial functions of fj and t)kw- 

We numerically solved the dilaton-gravity equations (jH)-((0J) using a standard Runge-Kutta 
routine. We adopted initial conditions around Eq ~ 1000, tq ~ 1, which comes from the 
requirement of T-duality. The shifted dilaton is chosen to satisfy the condition e* < 1 to ensure 
that the string coupling constant is initially small and hence that perturbation theory and the 
ideal gas approximation are trustable. The initial condition for A is somewhat arbitrary and we 
have carried out simulations for a wide variety of different initial values of A. Notice that ipQ 
is fixed by the constraint equation Q and that the negative solution is taken to remain in the 
perturbative regime of small string coupling constant. 

For initial conditions r$ G [0.8, 1.2] and Eq ~ 1000, the temperature is very close to the 
Hagedorn temperature with a nearly constant value. This is clear from eq. (|2H) . since the last 
term in (|21|) is vanishingly small relative to the first term, due to the exponential suppression 
given by e~^ E . Similarly we have P ~ for the above initial conditions from eq. I|22|) . Therefore 
the system is effectively described by a pressureless dust as shown in fig. ^ In this case one has 
A ~ Ae^ from eq. @, with A an integration constant. Subtracting eq. (jfjj) from eq. (jlj), we find 
a simple relation, (e -- ^) = E/2. Taking note that E is nearly constant (E ~ Eq), the analytic 
solutions of eq. (| H )-([B j) in the Hagedorn regime may be written as 



En o B 
-V + Bt + - 
4 



dA 2 



A 



Ao + ^log 



Eq 

(E t + 2B 



2VdA)(B + VdA) 



(E t + 2B + 2VdA)(B - \fdA) 



(23) 
(24) 



A and B are integration constants depending on the initial values for A, tp and ip. In particular 

A = \Qe-^°, B = -ii e-^°, (25) 

and d is the number of dimensions (we are now considering the case with d = D = 9). Notice 
that due to eq. (@J), tpQ, and thus B, can not be taken to be vanishing. 
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Figure 1: The evolution of r and P for the initial conditions Ao = 10~ 3 , Ao = 0.2, ipQ = —5 and 
Eq = 10 3 in the Hagedorn regime, ipo is determined by the constraint equation Ijl]). We plot the 
evolution of r both using the analytic approximation ()24l) and the full numerical result, which 
show very good agreement each other. 



In fig.^we plot the evolution of r that corresponds to the analytic solution (J24[). together with 
the full numerical results. They show very good agreement, which implies that the Hagedorn 
regime is well described by a state with a constant energy and negligible pressure. This actually 
ensures the validity of the analytic estimation in ref. ^21 discussed briefly in its Appendix. 

As long as Ao is positive (negative), the radius grows (decreases) towards the asymptotic 
value 



B + VdA 1/W 



B-VdA 



(26) 



with r getting smaller with time (see fig. ^) . We have checked this for values of ro very close to 
1, up to ro = (1 + 1 x 10~ 15 ), and found no substantial changes in behaviour. If one chooses 
exactly ro = 1, eq. (|15[) should be replaced by another similar relation, since now the two poles 
of order 18 in the /3-plane approach each other to a single pole of order 36. The above results 
apply also in this case: the pressure is almost zero and the evolution of the system is very slow 
in time. 

For initial values of E in the range Eq £ [500, 5000], the dynamics of the system is practically 
the same as explained above. For initial energies Eq > 1000 the scale factor is essentially 
constant in time. On the other hand the system is typically unstable for Eq <C 1000 and not 
thermodynamically meaningful for such low values of Eq. 
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3.2 Conservation laws 



We can also address the question of what happens when conservation laws are taken into account. 
In the case of the toroidal compactification we consider, the conserved quantities are taken to be 
the total winding number N and the KK momenta Mj in each compact dimension (i = 1, . . . , 9). 
This is performed by introducing a chemical potential for each conserved charge |3Uj . 

A crucial difference with respect to the previous case is that the leading singularity of the 
partition function, namely the Hagedorn temperature, is enough to study the thermal behaviour 
of the system. In particular, the leading singularity now depends on the compactification radius 
induced by the conservation laws. The entropy of such an ideal gas can be easily computed, 
yielding 

S(E, r) ~ p H E - ^ + N " r2 ) ~ 9 ln E ■ ( 27 ) 

We see that in eq. (|27(l there are two suppression terms with respect to the leading term (3hE in 
eq. ()15[). This is expected, because eq. (fTa|) counts states with all values of charges whereas (f2*7|) 
counts a smaller set of states in which the value of the charge is fixed. Moreover, the number of 
states decreases as Mj or N are increased, since less energy is available for the oscillators. The 
temperature and pressure obtained from the entropy (|27|) are: 

9 2 v / 2#^V r / 

The pressure vanishes if one imposes vanishing winding and KK charge, Mj = N = 0. The 
energy E evolves very slowly for Eq of order 1000 in which case the radius r asymptotically 
approaches a constant value after some growth from ro — 1, thereby showing similar behaviour 
to fig. E When Mi,N ^ 0, due to presence of the E factor in the denominator of eq. (|29j). the 
system evolves similarly to the case Mj = N = 0, as long as the summation terms in eqs. © 
and iJBJ) are unimportant. The dynamics changes if Mj and N are of order 10 5 , since the scale 
factor can have respectively a significant expansion or decreasing rate. As expected, winding 
modes prevent expansion, whereas KK modes, as standard matter, favorite it. 

As a last remark, notice that the string gas defined by eqs. (|2*S)l and ((2*9")) has a negative 
specific heat. Systems with negative specific heat are unstable in non-compact spaces, but 
actually can be in thermal equilibrium in a finite space. Along the lines of |31l I32| . we have 
evaluated the critical volume V c under which the system is in equilibrium and found that this 
is actually the case for V c S> 1, implying that the system can be actually in equilibrium with 
radiation. It should also be emphasized that this system is trustable as long as string corrections 
are completely neglected. When string interactions are included, the system is most likely to 
undergo a phase transition |33j whose details are so far unknown. 

4 Almost-radiation regime 

Assuming that some dimensions (d) start to expand while the remaining (9 — d) dimensions 
remain small by means of some mechanism, the system will eventually reach a temperature 
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below the Hagedorn regime where the dynamics is mostly governed by massless states, i.e. 
radiation 5 . In this case, it is important to see the conditions under which the large dimensions 
continue to expand. At the same time, it is important to study whether the small dimensions 
remain small or if they also enter an expanding phase. Different from the Hagedorn regime 
discussed in the section 3, the microcanonical ensemble agrees with the canonical one for the 
range of temperatures and energy densities involved now. In the following we shall use the latter 
ensemble, which is more convenient for practical purposes. 

Assuming again thermodynamical equilibrium and adiabatic evolution, we shall solve the 
dilaton-gravity equations P ]) -(|12j l . with E and P obtained from the free energy of a string gas 
at temperatures below Tjj- The entropy S is conserved with time under the assumption of 
adiabatic evolution: 

dt dt\ dp J v ; 

Eq. l|30j) is solved by letting [3 and the scale factors Aj = log Ri, be slowly varying functions of 
time \p — > P(t), Aj — ► Aj(i)]. In this way one can derive a differential equation whose solution 
gives /3 = /3(Aj) with S =constant. We denote the radii of the large d dimensions, taken all 
equal, by R = e^, whereas the radii of the (9 — d) small dimensions, again all equal, by r = e v . 



4.1 Pure radiation 

As a first step, let us consider the case of pure radiation (see the Appendix for the free energy 
in the context of the canonical ensemble). The energy and pressure are easily evaluated from 
eq. : 

e ( 1 = f%1 + /3^P = ^(o) 2 r ( *±1) (4^c(d+ i)(i - 2-( rf + 1 ))/r d - 1 , (31) 



d/3 2vr 

whereas the pressure P~^l for the d spatial dimensions is given by 



p( d ) _ 1 C 11 'rad _ _rpi d ) _ ^rad (on\ 

rad dd()nR) rad d ' 



which corresponds to the equation of state for radiation in d spatial dimensions. Eq. (|31[) is 
nothing but the d-dimensional generalization of the Stefan-Boltzmann law in presence of D(0) 2 /2 
bosonic and fermionic degrees of freedom. Since F^ d does not depend on r, the pressure along 
the small dimensions vanish: 



p(9-d) _ 1 d F rad . n ,•>•-)> 

Kad ~ 9 _ dd{ l Qr) - U - , • > ' ,, 



From the adiabatic equation (|3(jp. we easily get the following relation 



P = Po^~, (34) 



5 Notice that the system is already in an almost radiation regime for (3 > 11 and it is essentially governed by 
pure radiation only for j3 > 14, in string units. 
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relating the temperature and scale factor in a radiation-dominated universe, with /3q and Rq 
being initial conditions satisfying (3{Rq) = (3q. In this case the dilaton-gravity equations ([§])- 
(O read 

^ = ^ + l {9 ~ ^ + \^ ' (35) 
ii = ^ A +I e *pW, (36) 

P = -0^ ) (37) 

together with the constraint equation 

V> 2 = ^4ad + d A 2 + (9 - d)i? . (38) 

Eq. (|3*7|) is integrated to give 

v = 2> oe ^° , (39) 

where z>o are 'i/'o are the initial values of v are ip. We see from eq. (|39|) that when z>o is positive 
(negative) the expansion rate for the small dimensions is always positive (negative). In order 
to avoid unbounded growth of the dilaton towards the strongly coupled regime (e* > 1), it is 
natural to consider the case with negative ip. In this case the absolute value of v decreases with 
time. 

In the absence of the pressure P ^ X in eq. (|36|) , the evolution of the large dimensions is similar 
to that of the small ones. In the case ipo < and fiQ > 0, we have ft < for Pzfl = from 
eq. (|36|) . This corresponds to the universe with expanding large dimensions with a decreasing 
Hubble rate. Numerically we found that the evolution of the system in this case is trivial, namely 
the large dimensions soon approach a nearly constant value with very small /i. 

When the pressure P^ d is taken into account, this works as a positive source term in eq. (|36|) . 
Therefore it is possible to make the r.h.s. of eq. (|36|) positive even when ip < and /i > 0. We 
have made numerical simulations with initial conditions Rq 3> vq ~ 1, ipo < 0, and several 
different values of fiQ and uq. As long as /io is positive, the large dimensions expand in the 
presence of the pressure due to radiation. The contribution of the pressure term in eq. ()36|) 
inhibits the rapid decrease of /i, thereby leading to different evolution of R compared to the 
case with P„S = 0. The expansion rate v for the small dimensions is exponentially suppressed 
with the decrease of ip [see eq. ()39|) ]. Therefore unless the initial value of \i>\ is much larger than 
unity, the radius r can stay small around r ~ 1. 

We have numerically succeeded to obtain ideal solutions with growing R and small, roughly 
constant, r satisfying r <^ R- One typical evolution is plotted in fig. El These solutions can be 
achieved by choosing initial values with /io > l^oli Rq ^> r o ~ 1 and ipo < 0. When /io "C l^ol 
holds initially, it is difficult to keep the small dimensions small relative to the large ones. If 
/io < 0, we have /io > from eq. (|36]> . This leads to the growth of the expansion rate /i. Since 
ji continues to be negative by the time it crosses zero, the large dimensions contract during this 
stage. After /i changes sign, R begins to grow. This implies that bouncing solutions may be 
obtained if /io < 0. We have numerically found that this is actually the case, see fig. |3 During 
the contracting phase, the temperature increases according to eq. I|34|) . The temperature is 
maximum at the bounce where R is minimum. In the context of Pre-Big-Bang ^ or Ekpyrotic 
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20 40 60 80 100 

t 



Figure 2: The evolution of R, r and (5 for the pure radiation case with d = 3. We choose the 
initial conditions /io = 1.0, //o = 4.0, z>q = —0.01, uq = 0.0, Vo = —16 and /3q = 15. 




Figure 3: The evolution of R, r and (5 for the pure radiation case with d = 3. We choose the 
initial conditions /io = —0.04, = 4.0, uq = —0.01, vq = 0.0, = —16 and /3q = 15. 
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Cosmologies jS] , nonsingular bouncing solutions are difficult to construct unless loop or derivative 
corrections are added to the tree-level action |34| I35| . It is quite interesting to be able to obtain 
bouncing solutions only by including radiation in dilaton- gravity equations. 

It is worth investigating the asymptotic behavior of the dynamical system of eqs. (|35 [l -(|37 [) . 
along the lines of |12j . Introducing new parameters, fi = £, i> = r], tp = f, and using the fact 
that the pressure in the large dimensions is connected to the energy Ef} d through eq. 
eqs. (|33|) - (|3*%|) become: 



/ 

i 

V 



e + 



i n< 



1 

fv- 



e + 



1 \rf 



2d 



(40) 

(41) 
(42) 



Since we are considering the case with decreasing if), rj asymptotically approaches zero from 
eq. (|39|) . i.e. r] = is an attractive solution. In this case, the analysis is closely related with the 
one outlined in the Appendix of ^2]. In particular one finds that the line described by 



ff£ = -d, 77 = 0, 



(43) 



is an attractor. From eq. (j5J) the time-derivative of the dilaton is given as 2$ = f + d£ + (9 — d)rj. 
Therefore we have $ = for the attractor (|43jl. again in complete analogy to the case of the 
single scale factor |12| . Substituting eq. (|43|) for eq. (|41|) and integrating this relation, one finds 



(d+l)t' 



R oc . 



(44) 



This indicates that the late time evolution for the large dimensions can be described by that of 
the standard radiation dominant phase in FRW cosmology for d = 3, even in the presence of 
the small dimensions. The key point is that the pressure in the small dimensions vanishes for 
the massless case, thereby leading to jj = as an attractor. Notice that a cosmological solution 
of this kind has been also obtained in standard General Relativity and in a purely Kaluza-Klein 
extra-dimensional scenario by |26j . 

We also analyzed the evolution of the system by varying the value of d, and found that the 
situation is not basically changed compared to the d = 3 case discussed above. As long as the 
initial conditions satisfy fio > |z>o|, Ro 3> r$ ~ 1 and ipQ < 0, the large dimensions continue to 
expand due to the presence of radiation while the small dimensions are kept to be small (r <C -R). 
In this case the large dimensions asymptotically approach the radiation-dominant FRW solution 
(|44|). When the initial value of /i is negative, we found that it is possible to have a bouncing 
cosmological solution that approaches the expanding FRW universe given by (|44jl. 



4.2 Inclusion of matter 

Although massless states, pure radiation, dominate the thermodynamical ensemble in this phase, 
this dominant contribution has a trivial dependence on the small dimensions r. In particular, 
iyfj does not depend on r and the pressure along the small dimensions trivially vanishes. It is 
then important to see if and how matter terms can alter this behaviour. For this purpose, we 
study the leading terms that have an explicit dependence on r in the infinite sums appearing in 
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eq. (|54|). We have numerically estimated that it is enough to consider the first KK and winding 
modes along a small direction, i.e. the terms with {iV = N = 0, mi = (1,0,... , 0),nj = 0} (as 
well as rrii and n% exchanged) in eq. (|54|). plus the remaining 8 — d inequivalent permutations. 
The energy E^ at and pressures -P^fat an d Pmat^ along the large and small dimensions associated 
with these states are simply evaluated starting from the general expression (|54|) . The equation 
of state for these leading order terms are: 



P, 



mat 



(9-d) 
mat 



r (d+i)/2 2 ( d+1 )/ 2 yTj r (d+i)/2~ n (d+i)/2 y- 

'r^ 2 ^K (d+1)/2 (Pr) + r^(3rK[ d+1)/2 ((3r) ' 
^^(/ 3 ) (d) [^^^ + l)/2(7)+^ d+1)/2 ^ + l)/2(^)^ 



+ 



(45) 
(46) 
(47) 



where K n are modified Bessel functions, the prime denotes derivative with respect to r, and 



C(J3) 



(d) 



(d+i)/2 

0/ 



2d) 



7T 



D(0) 5 



(48) 



Here L>(0) 2 = 256 is a string degeneracy factor (see the Appendix). Note that we only consider 
the p = 1 term in eq. (|54|). implying the approximation of the bosonic/fermionic statistics with 
the Maxwell-Boltzmann distribution. The pressure -P^at along the large dimensions is always 
positive, which aids expansion of the universe in addition to the pressure P~gl from the massless 
states. The first and second terms in square brackets in eqs. (|46[) and (|47|) come from the KK and 
winding mode, respectively. The above equations (|45 |) -(|47 |) are all manifestly invariant under 
the duality symmetry Q acting on the small dimensions, r — ► 1/r. Notice that the winding 
modes give rise to a standard positive pressure along the large dimensions (second term in (|46[)) 
but negative along the small ones (second term in (|47j0 . 



We numerically solved the dilaton-gravity equations P )l -(|12 |) . with E = E, 



J 1 ~T~ -t 77 



rad ± mat 

condition (|3fl|) 



and P^-^ 



p (9-d) p (9-d) 
rad ' " 



(d) 
rad 



+ E (d) p{d) 
~ 'mat i ± 



The pressure P^at^ f° r the small dimensions vanishes at the self-dual critical 
radius r = 1. Therefore it is expected that the effect of the massive states for the small 
dimensions is weak around r ~ 1. In fact we have numerically found that this is the case. 
As seen from the case (b) in fig. the evolution of the small dimensions is hardly altered by 
including the massive mode for the initial value of r very close to unity. From eq. (|47|) one notes 



mat 



by carefully taking into account the adiabaticity 



that Pjnat < for < r < 1 and Pm at a> > for r > 1 (the asymptotic values are Pmat U ' ~^ 
for r — > and r — > oo). This indicates that the pressure of the massive state makes the small 
dimensions contract for < r < 1 while its effect tends to expand the small dimensions for 
r > 1. 

The effect of the massive states emerges by choosing the initial values of tq that are slightly 
smaller or larger than unity. For < tq < 1 with z>o > 0, the small dimensions can be larger 
than r = 1 for large initial values of v. In this case the small dimensions continue to grow after 
they cross r = 1. When z>o is not large (z>o <C 1), the massive effect can lead to the contraction 



,(9-d) 



(9-d) 



14 




Figure 4: The evolution of the small dimensions for d = 4 when the massive states are taken into 
account (solid curve). The dotted curves correspond to the case where the massive states are 
neglected (only massless states). We choose the initial conditions /to = 1-0, (jlq = 4.0, ipo = —16, 
Po = 12 with (a) z> = -0.1, u = 1.0, (b) z> = -0.1, u = 0.05, (c) z> = 0.1, u = -1.0. 

of the small dimensions due to the negative pressure for r < 1. As found from the case (c) in 
fig. the small dimensions always increase in the massless case, whereas the small dimensions 
begin to contract if the massive effect is included. Therefore we can keep the (9 — d) dimensions 
small (0 < r < 1) for these initial conditions. 

We have also made numerical simulations for tq > 1. When (iq > and z>o > 0, both large 
and small dimensions expand in the presence of positive pressures. If vq is largely negative, 
the small dimensions contract by passing through r = 1. Meanwhile, if <C 1, the small 
dimensions can exhibit bouncing with r > 1, instead of crossing r = 1 [see the case (a) in 
fig. BJ. This means that the radius r can grow in the presence of the massive states. Since the 
small dimensions continue to expand after the bounce, this is not an ideal case where the small 
dimensions stay small. Nevertheless the small dimensions can be made small compared to the 
large dimensions as long as £iq > Vq. 

When the massive states are taken into account, this gives rise to an extra source term for 
the energy E in eq. ©• Then |?/>o| gets larger compared to the pure massless case. Typically 
this leads to the suppression of the r.h.s. of eq. (|10j) via the term, thereby yielding the smaller 
expansion rate /t for the large dimensions. The large dimensions in the massive case grow slower 
relative to the massless case, as seen in fig. 03 In addition, the massive effect suppresses the 
growth of j3, i.e. the temperature decreases faster in the massless case. 

As expected, the massive terms get smaller as the initial value of (3 is increased. For example, 
in the cases shown in fig. 21 and EJ the system is effectively described by the massless states for 
Po > 15. We also analyzed the behaviour of the system by varying the value d with 1 < d < 8. 
We found that the numerical results are quite similar to the case explained above (d = 4). As 
long as the conditions, /to > |^o| and Ro > r o ~ 1; are satisfied, the small dimensions are kept 
small, while the large ones expand as required in the presence of radiation and string matter. 
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-20 20 40 60 80 100 

t 

Figure 5: The evolution of the large dimensions and (3 that corresponds to the case (c) in fig.EJ 
The dotted curves correspond to the case where the massive states are neglected (massless states 
only). 

We have also considered the case with the conservation of KK and/or winding modes. Under 
this circumstance, the partition function must be taken summing only over the configurations 
that respect the conservation laws. As before, this is done by introducing a chemical poten- 
tial corresponding to each conserved quantity in the partition function. It turns out that no 
significant changes are found compared to the case where no conservation laws are imposed. 

5 Conclusions 

We have studied string-gas cosmology at finite temperature in a toroidal universe. We make use 
of the dilation-gravity equations of motion, satisfying the R—>l/R duality symmetry, to study 
the evolution of the system, which is assumed adiabatic. Our set up is as follows: the universe, 
initially homogeneous, isotropic and in thermal equilibrium, is filled with an ideal gas of closed 
strings. The 9 spatial dimensions, compactified on a 9-torus, evolve adiabatically starting from 
a Hagedorn regime. 

In order to find an acceptable equation of state that describes such a system in string theory 
at finite temperature, we evaluate the energy and pressure in a microcanonical approach. Close 
to the Hagedorn regime, the scale factor R exhibits a slow time evolution around R ~ 1, as 
shown in fig. ^ I n this case the dynamics of the system is effectively described by a nearly 
constant energy and negligible pressure. We found that the analytic solution in the Hagedorn 
regime shows very good agreement with the full numerical result. No substantial changes are 
observed even when the conservation of Kaluza-Klein (KK) and winding modes is imposed, as 
long as the conserved charges are of order one. 

We have also investigated a "low" temperature regime, in which the equation of state is 
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derived in a canonical context. We first considered the dynamics of 3 "large" and 6 "small" 
compact dimensions in the presence of a pure gas of radiation (given by the massless states). 
It turns out, as expected, that there exist interesting cosmological solutions where the large 
dimensions continue to expand while the small dimensions remain nearly constant and small 
relative to the large ones (see fig. |2J). The attractor solutions for the large dimensions can 
be described by the evolution of the standard radiation dominant phase in FRW cosmology 
whereas the small dimensions always asymptotically approach to a constant value. We also 
found bouncing solutions for the large dimensions if their Hubble rates are negative initially (see 

fig.H. 

We then analysed the case where the massive string states are taken into account in addition 
to the pure radiation. The presence of the massive states typically leads to a slower expansion 
of the "large" dimensions relative to the massless case (see fig. EJ). Meanwhile the behaviour of 
the small dimensions strongly depends on the initial conditions for r and r, resulting in either 
expansion or contraction of the small dimensions (see fig.|3J). The radius r can be kept small as 
long as r is initially close to unity, since the pressure vanishes at the duality symmetric radius 
(r = 1). The vanishing of the pressure at r = 1 is a purely stringy effect, since it is due to 
winding modes, whose negative contribution compensates that of KK states. The important 
point is that, even in the presence of the massive state, there exist a wide range of the initial 
condition space for which the small dimensions are stabilized around the self-dual radius and 
are kept small relative to the large ones. These behaviours are found to be insensitive to the 
number of large dimensions, d. We also considered the case for the conservation of KK and 
winding modes and found no substantial change compared to the case without imposing the 
conservation laws. 

In this work we did not investigate the later stage of the cosmological evolution. From the 
view point of cosmology, it is important to have an inflationary stage in order to solve the 
major cosmological problems (such as the flatness problem) as well as to generate the seeds for 
large-scale structure. One way to address this problem is to assume that the dilaton (or moduli) 
acquires mass, which may lead to inflation at later stages. Recently, for example, an interesting 
proposal was made by Parry and Steer [201, wno showed that inflation can occur on a moving 
3-brane due to the nonminimal coupling of the dilaton to the brane matter. It was also shown in 
ref. that brane inflation could result from the negative pressure of winding modes. Although 
it is not clear at present whether these proposals are directly applicable to our string-gas work, 
it is certainly of interest to extend the present analysis in that direction. 
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APPENDIX : The one-loop string partition function 



where 



A thermal canonical partition function at temperature T can be computed in the Euclidean 
approach by compactifying the time direction on a circle of length (5 = 1/T. We are then 
led to consider Type IIA/B string theory with all ten space-time directions compactified. For 
simplicity, we take the tori to be all rectangular (i.e. simple products of circles) and compute 
the free energy for the case in which (9 — d) spatial dimensions are compactified on small radii all 
equal to a common value denoted r, whereas the remaining d directions are taken very large, and 
all equal to R. The free energy, after having unfolded the fundamental region of the world-sheet 
torus, takes the following form (see e.g. 38 ): 

FW(/3) = --4= f\d n r ^^[Hr ]T )r d ±e-^\M 2 \\r), (49) 



2 T 2 p=1 

+oo , 

A(r;r)= ]T ql(f+^- W -nr )^ (5Q) 
m,n=— oo 

represents the contributions of the whole KK and winding modes along the small dimensions, r is 
the so-called modular parameter of the world-sheet torus and q = exp(2i7rr). In deriving eg. (1491) 
the winding modes along the large dimensions have been completely neglected and the sum over 
the KK modes has been approximated by an integral over continuous momenta. Note that p 
in ea.(|49|) runs only over positive odd numbers and correspond to taking the correct quantum 
statistic for bosons and fermions. Taking only the term p = 1 in the above sum corresponds 
to replacing the quantum bosonic/fermionic distribution with the classical Maxwell-Boltzmann 
distribution. Vd is the volume of the large dimensions in (47r 2 ) units 6 : 

V d = - * (2nR) d = R d , (51) 



( 47r 2)d/2 

The M2 factor in eq. (|49|) encodes the contribution to the free energy of the whole tower of 
massive strings, and can be expanded in powers of q: 

^) = ^=±D(N)q». (52) 

Here #2 and r\ are modular functions on the torus (see e.g. jSH] for an explicit expression) and 
D(N) is the degeneracy factor at level N (D(0) = 16, for example). The value of iV corresponds 
to each string mass level. The r\ and T2 integrals in eq. Q49j) can be easily performed. It is 
convenient to consider the term with N = N = rrii = rii = in eq. ()49|) separately from the 
remaining ones. This is the contribution of the purely massless states, which we will henceforth 
denote as F ra( i (where "rad" stands for radiation). We get 

4al = "i^(°) 2r (^) ^) {d+1)/2 «d+ 1)(1 - 2^+ 1 ))/r d - 1 , (53) 



3 Recall that we are setting a' = 1. 
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where £(x) is the Riemann zeta-function. The remaining "matter" terms give 



51(d) 

mat 



— E E E , " 2kM ^' i+ " > 1 



.LLhr ^n i +iV-JV,0^(^)^(^)^(rf+l)/2(^M), (54) 



mi,rii p=i jv, JV 

where if(d+i)/2 are modified Bessel functions, and 



M = M{m h ni,n,N,N) 



\ 




+ nfr 2 + 2(JV + N) . 



The total free energy is given by 



F (d) F (d) „ 



rad 



mat 



(55) 



(56) 



The infinite sum over N and N is not always convergent. In fact, the degeneracy factors D(N), 
for large values of N, have a leading exponential behavior D(N) ~ exp(27ry / 2iV). On the other 
hand, for large values of its argument, the modified Bessel function K n (z) admits an asymptotic 
expansion whose leading term is ~ exp(— z). Hence, the sum over N and N in eq. (|54(1 converges 
only for (3 > /3h = 2%y2. The temperature Tjj = 1/Ph is the Hagedorn temperature. 

As long as we deal with a range of energies where eq. (|49[) converges and no large energy 
fluctuations are present, we can work with the canonical ensemble. On the other hand, for high 
energy densities a microcanonical description has to be used. In this case the energy density of 
states is governed by the analytic structure of the canonical partition function Z = exp(—f3F), 
in the complex /3-plane. Taking into account the leading singularities of (|49jl. we can parametrize 
the partition function as 



Z(P,R) 



P-Ph\P-P. 



vk 



7K 
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Vw 



P-Pw 



IS 



(57) 



where t)k and rjw are defined by eq. (fTB|) with fix = Ph — i]K and Pw = Ph — f]w- R) is 
an entire function in p. The microcanonical energy distribution function £l(E), is then given by 



f3 H +ioo f]Q 



pH—ioo 



2ni 



i=H,KW 



(58) 



where Ch, Ck and Cw are the three contours encircling respectively the poles in pjj, Pk and 
Pw in the complex /3-plane. The entropy and the rest of the thermodynamical quantities easily 
follow from eq. (|58|) . 



References 

[1] I. Antoniadis, Phys. Lett. B 246 (1990) 377; 

I. Antoniadis, C. Munoz and M. Quiros, Nucl. Phys. B 397 (1993) 515 
|arXiv:hep-ph/9211309| . 

[2] R. Brande nberger and P. M. Ho , Phys. Rev. D 66 (2002) 023517 [AAPPS Bull. 12N1 

(2002) 10] |arXiv:hep-th/0203119| ; 

S. Alexander, R. Brandenberger and J. Magueijo, arXiv:hep-th/0108190 ; 
S. Alexander and J. Magueijo, |arXiv:hep-t h/0104093 



19 



[3] J. Maldacena, arXiv:astro-ph/0210603 



[4] G. Veneziano, Phys. Lett. B 265, 287 (1991); 

K. A. Meissner and G. Veneziano, Phys. Lett. B 267, 33 (1991). 

[5] M. Gasperini and G. Veneziano, Astropart. Phys. 1 (1993) 317 |arXiv:hep-th/9211021| . 

[6] J. E. Lidsey, D. W ands and E. J. Copeland, Phys. Rept. 337, 343 (2000) 
|arXiv:hep-t h/990906l"l . 

[7] M. Gasperini and G. Veneziano, arXiv:hep-th/0207130 . 

J. Khoury, B. A. Ovru t, P. J. Steinhardt and N. Turok, Phys. Rev. D 64, 123522 (2001) 
|arXiv:hep-th/0103239| ; Phys. Rev. D 66, 046005 (2002) arXiv:hep-th/0109050| . 



[9] P. J. Steinhardt and N. Turok, arXiv:hep-th/0111030 ; Phys. Rev. D 65, 126003 (2002) 
|arXiv:hep-th/0111098] . 

[10] H. Liu, G. Moore and N. Seiberg, JHEP 0210, 031 (2002) |arXiv:hep-th7 0206182 ; 
S. Elitzur, A. Giveon, D. Kutasov and E. Rabinovici, JHEP 0206, 017 (2002) 
|arXiv:hep-th/0204189] . 

[11] R. H. Brandenberger and C. Vafa, Nucl. Phys. B 316 (1989) 391. 

[12] A. A. Tseytlin and C. Vafa, Nucl. Phys. B 372 (1992) 443 |arXiv:hep-th/9109048| . 

[13] A. A. Tseytlin, Class. Quant. Grav. 9 (1992) 979 |arXiv:hep-th/9112004| . 

[14] M. A. Osorio and M. A. Vazquez-Mozo, Mod. Ph ys. Lett. A 8 (1993) 3111 
|arXiv:hep-th/9305137| ; Mod. Phys. Lett. A 8 (1993) 3215 |arXiv:hep-th/9305138| . 

[15] G. B. Cleaver and P. J. Rosenthal, Nucl. Phys. B 457 (1995) 621 |arXiv:hep-th/9402088| . 

[16] M. Sakellariadou, Nucl. Phys. B 468 (1996) 319 |arXiv:hep-th/9511075| . 

[17] S. Alexander, R. H. B randenberger and D. Easson, Phys. Rev. D 62 (2000) 103509 
|arXiv:hep-th/0005212] . 

[18] R. Brandenb erger, D. A. Easson and D. Kimberly, Nucl. Phys. B 623 (2002) 421 
|arXiv:hep^th/0109165l. 



[19] D. A. Easson, arXiv:hep-th/01 10225 



[20] R. Easther, B. R. G reene and M. G. Jackson, Phys. Rev. D 66 (2002) 023502 
|arXiv:hep-t h /0204099] . 



[21] S. Watson and R. H. Brandenberger, arXiv:hep-th/0207168 



[22] R. Easther, B. R. Greene, M. G. Jackson and D. Kabat, arXiv:hep-th/0211124 
[23] Y. Leblanc, Phys. Rev. D 38 (1988) 3087. 



[24] T. Damour and A. M. Polyakov, Gen. Rel. Grav. 26, 1171 (1994) arXiv:gr-qc/9411069 . 



20 



[25] T. Damour, F. Piazza and G. Veneziano, Phys. Rev. D 66, 046007 (2002) 
|arXiv:hep-th/0205111| ; Phys. Rev. Lett. 89, 081601 (2002) |arXiv:gr-qc/0204094| . 

[26] S. Randjbar-Daemi, A. Salam and J. Strathdee, Phys. Lett. B 135 (1984) 388. 

[27] L. Alvarez-Gaume and M. A. Vazquez-Mozo, arXiv:hep-th/9212006 , section 4.3 and refer- 
ences therein. 

[28] R. Hagedorn, Suppl. Nuovo Cimento 3 (1965) 147. 
[29] N. Deo, S. Jain and C. I. Tan, Phys. Lett. B 220 (1989) 125. 
[30] N. Deo, S. Jain and C. I. Tan, Phys. Rev. D 40 (1989) 2626. 
[31] S. W. Hawking, Phys. Rev. D 13 (1976) 191. 

[32] M. Axenides, S. D. Ellis and C. Kounnas, Phys. Rev. D 37 (1988) 2964; 
M. J. Bowick and L. C. Wijewardhana, Phys. Rev. Lett. 54 (1985) 2485. 

[33] J. J. Atick and E. Witten, Nucl. Phys. B 310 (1988) 291. 

[34] M. Gasperini, M. Maggiore and G. Veneziano, Nucl. Phys. B 494, 315 (1997) 

|arXiv:hep-th/9611039| ; 

R. Brustein and R. Madden, Phys. Rev. D 57, 712 (1998) |arXiv:hep-th/9708046| . 

[35] S. Tsujikawa, R. Bra ndenberger and F. Finelli, Phys. Rev. D 66, 083513 (2002) 
|arXiv:hep-th/0207228] . 

[36] M. F. Parry and D. A. Steer, JHEP 0202, 032 (2002) |arXiv:hep-ph/0109207| . 

[37] S. Abel, K. Freese and I. I. Kogan, JHEP 0101, 039 (2001) [archiv; hep^th/0005028| ; 
|arXiv:hep-t h /02053 17} 

[38] E. Alvarez and M. A. Osorio, Phys. Rev. D 36, 1175 (1987); 
K. H. O'Brien and C. I. Tan, Phys. Rev. D 36 (1987) 1184. 

[39] J. Polchinski, String Theory, Volume 1, Cambridge, page. 214-216. 



21 



